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^ ! Abstract A brief introduction to the String Fragmentation Model and its consequences 
^ I is presented. We discuss the fragmentation of a general multi-gluon string and show that 
^ ' it can be formulated as the production of a set of "plaquettes" between the hadronic 
^ ■ curve and the directrix. We also discuss certain interesting scaling properties of the 
partonic states obtained from standard parton shower algorithms like those in Pythia 
^ ! and Ariadne, which are communicated to the final state hadrons obtained through our 
hadronisation procedure. 

1 The Area Law and its Consequences 

In the Lund String Model the massless relativistic string is used to model the QCD field 
^ ! between coloured objects, [jl|, 0, 0. The string has a constant tension k, which gives rise 
^1 to a linear potential. The endpoints of the string are identified as the quark and the 
anti-quark. The endpoints will stretch the string when they move apart and therefore 
lose energy-momentum. Because of this they will sooner or later have to reverse their 
direction of motion and again gain energy-momentum from the string. The motion of a 
simple string system with a massless quark and anti-quark at the endpoints will therefore 
be the so-called yoyo- mode, ||, ||. We will now describe this string state in more detail. 

In order to be more specific, let us consider the e"'"e~-annihilation process in Fig.[^. 
We assume that the quark and the anti-quark are produced at a single space-time point 
and that they start to move apart. This means that there will be one original string 
spanned between the quark and the anti-quark, cf. Fig.|l]. Gluons are interpreted as 
internal excitations of the string and they will complicate the string motion dramatically. 
In this section we will only consider string-states without gluons, and we will come back 
to multi-gluon strings in the next section. The final state hadrons observed in high energy 
e"'"e~-collisions stem from the breakup of the force field. New quark anti-quark pairs are 
produced from the field. This means that the original string will be split into smaller 

^Thc article is compiled by S. Mohanty and F. Soderberg from the transparencies of the talk presented 
by the Late Prof. Bo Andersson. 
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Figure 2: Two adjacent vertices and a hadron with energy-momentum pj produced in 
between. 

string pieces because the quarks and anti-quarks are identified as string endpoints in the 
string Model We are now going to describe this break-up process in more detail. 

A typical Lund String break-up process can be found in Fig.|l]. A set of new gg-pairs 
are produced at certain space-time points (called vertices). One can note the following 
things in this figure: 

• There is no field between a gg-pair produced at a single vertex. A string force field 
is always confining because it has a fixed energy per unit length and the force field 
vanishes at the end-point charges. 

• The hadrons must have positive momenta which implies, cf. Fig.|l], that the vertices 
must be space-like separated. This means that no invariant time-ordering is possible. 
There is no principle "first" vertex, they are all "equal". In the Lund Model, it will 
always be the slowest particles that are produced first in any frame. 

1 

A convenient ordering prescription is rank-ordering. Two hadrons that share the same 
vertex will be defined to have adjacent rank. If one starts to count from the original 
quark end, this means that the hadron which contains this quark and an anti-quark from 
a break-up vertex will have rank 1, cf. Fig.[^. The break-up process can be described as a 
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Figure 3: The vertex V will subdivide the system into "right-movers" {pi, . . . ,pj) and 
"left-movers" {pj+i, . . . ,Pn). 



set of steps along the positive light-cone momentum, where each hadron takes a fraction 
Zj of the remaining light-cone momentum, cf. Fig.^. 

It is possible to calculate the probability for a string break-up inside the Lund String 
Model using the basic properties that we just have mentioned if one also includes a 
"saturation assumption", [|l], 0, |^. This additional assumption basically means that even 
if the energy of the original pair will increase without limit the distribution of the proper 
times of the production vertices will stay finite. The two main results are: 

• The probability for producing a hadron with mass m taking a fraction z of the 
remaining (positive) light-cone momentum is given by f{z)dz, where 

/(^) = N-il-zTexpi ) (1) 

f is called the Lund Symmetric Fragmentation Function (probability distribution). 
We note that small z values will be exponentially suppressed (with i), while large 
z values will be power suppressed. The final result is that the vertices typically will 
be along a hyperbola in space-time in the Lund Model. 

• The Area Law, which states that the (non-normalised) probability to produce a set 
of hadrons with momenta pi,p2, ■ ■ ■ ,Pn and total momentum Ptot is given by: 

n n 

dPniiPjh Ptot) = n Njd%5ip/ - m.^iY^p, - Ptot) exp{-bA) (2) 

, i.e. the phase-space times a negative exponential of an area. The area is the area 
spanned by the string state before it decays, cf. Fig.^ 

Up to this point, we have described the break-up process in space-time. We are now going 
to show that one can also describe the same process in momentum space. In order to do 
this, we first note that a vertex (V) subdivides the state into {pn, ■ ■ ■ ,Pj+i} © {Pj, ■ ■ ■ ,Pi}, 

V ' v ' 

"to the left" "to the right" 
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Figure 4: The production process described in terms of momentum transfers. 



cf. Fig.|^. From this figure we conclude that: 



(3) 



We can therefore describe the process according to the left-hand figure in Fig.^, if we 
define the momentum transfer q = {X~^ , —X~). This process can be repeated for any 
vertex. This means that String Fragmentation can be described as a ladder diagram like 
the one in the right-hand figure of Fig.|[ There is a duality relationship between the 
space-time picture and this last diagram, cf. the left-hand figure of Fig.|^. 
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Figure 5: One possible way to partition the area of a fragmenting string into contributions 
for each hadron. It shows the connection between the area in the area law and the area 
between the directrix and the hadronic curve. 



2 Fragmentation of Gluonic Strings 

Presence of internal gluonic excitations cause the string to trace complex curved surfaces 
in space-time. Approximating the partonic state with a finite number of gluons with 
transverse momenta above a certain scale is equivalent to approximating the curved sur- 
face of a string by a set of planar regions with the (scaled) parton energy momenta at 
the boundaries. The dynamics of such a string is most conveniently described in terms of 
the action principle which implies that the string world surface should be a minimal area 
surface. 
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This means that the details of the string world surface are completely determined 
by a specification of the boundary, i.e., the quark or anti-quark world lines. For an 
open string these two world lines differ by just a translation, and hence it is sufficient to 
specify only one of them. This boundary curve is called the directrix. The position of 
an arbitrary point X{t, a) on the surface of the string (parametrised by time 't' and the 
energy a between the point and the quark end point) is related to the directrix A{t) in 
the following way. 

X{t,a) = ^(Ait + a)+Ait-a)) (4) 

The function A{t) has the following periodicity and symmetry properties for a string 
which starts at a single point: 

Ait + 2E) = Ait) + 2P (5) 
A{-t) = -Ait) 

We see from Eq.(^ and Eq.(^ that the string surface is described by a left moving 
and a right moving wave, which bounce at the endpoints. It can also be seen that the 
internal excitations would move and affect the end points in colour order. The directrix 
of a string stretched between a certain set of partons is therefore obtained by laying out 
the energy-momenta of the quark, the gluons and the anti-quark in colour order. 

The minimal area nature of the string world surface also means that properties of 
hadrons created during string fragmentation are stable against small perturbations on 
the boundary, i.e. that the model is infrared stable. 

Since the directrix describes the string state, hadronisation of a string can be formu- 
lated as a process along the directrix. In ^ we presented one such formulation of string 
fragmentation, and we briefly mention some of the main ideas here. 

In m , we defined an iterative process along the directrix to describe string fragmenta- 
tion. We defined (abstract) dynamical variables called "vertex" vectors Xi. These variables 
would correspond to true space-time locations for production of quark anti-quark pairs 
if the string had no gluonic excitations. Each particle energy-momentum vector pi is 
constructed as a linear combination of a vertex vector and a certain section of the 
directrix fcj. The same pair of vectors Xi^i and ki also determines the next vertex vector 

1 z^x^ 

Pi = ZiXi_i + -(1 '—Y^)ki (6) 

2 mf 

X^={l- Zi)Xi-i + hi + ^^^)ki 

2 rrii 

Since with this construction, Xi^i + ki = pi + Xi, we associate a "plaquette" bounded by 
those four vectors with the particle indexed "i" . Each of the vectors Xj is shared between 
two plaquettes, and hence production of a set of hadrons from a given string or directrix 
means construction of a set of connected plaquettes along the directrix in this model. The 
initial value xq, for the series of vertex vectors, is taken to be the first segment of the 
directrix. The numbers Zi are related to the area of the plaquettes. 
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There is a one to one correspondence between the areas of the plaquettes constructed 
in this way and certain areas which can be associated with the production of individual 
hadrons in the string fragmentation model without gluonic excitations. Therefore, the 
sum of the areas of the plaquettes, which is also the area between the directrix and the 
"hadronic curve" (i.e. the curve obtained by placing the hadron energy-momenta one 
after the other in rank order. We called this curve the X-curve in [^), could be used 
as the area in the Lund Model area law. Therefore, area law on the whole event can 
be implemented by weighting the construction of each plaquette according to its area, 
which amounts to choosing the variables "Zj" in Eq. (P) according to the Lund symmetric 
fragmentation function. 

It was shown in ^ that in the limit of infinitesimal steps fcj along the directrix, the 
hadronic curve approaches a continuous curve in space-time which we called a P-curve. 
This could be thought of as the equivalent of the average hyperbolical fragmentation 
region in a simple gluon-less string fragmentation scheme in this model. There is a close 
connection between the P-curve and a certain A'-curve defined on the partonic state, or 
the directrix. The A'-curve and its length, called the A-measure which is related to the 
available phase space for hadron production, are easily calculated for a given directrix. 
It is apparent therefore that the average hadronic curve in this model is determined by 
perturbative calculations. 



3 Generalised Dipoles 

In this section we will discuss certain properties of the QCD bremsstrahlung mostly in 
connection with the dipole cascade model as implemented in the Monte Carlo routine 
Ariadne and relate them to the hadronisation procedure we have just mentioned. 

Gluon radiation from QCD colour dipoles can be compared to photon radiation from 
electromagnetic dipoles. The inclusive density for the emission of a gluon from a quark 
anti-quark dipole is given by 



dn = tteflf ■ -jY' ' '^y ' (PolarisationSum) (7) 

This is very similar to the result for photons. Energy momentum conservation implies 
that the transverse momentum and rapidity of the emitted particle obey the following 
constraint: 



k^coshiy) < (8) 

where ^/s is the center of mass energy of the dipole. This in turn, means that the 
total rapidity range available for the emission is given by 



Ay ^ In l^-^J (9) 

The emission of a photon leaves the electromagnetic current unchanged except for 
small recoil effects, whereas the emission of a gluon, which has colour charge, changes the 



6 




Figure 6: The triangular Lund phase-space diagrams described in the text. 

current for further gluon radiation. However, it was shown in [0, ^ that to a very good 
approximation, the density for the emission of two gluons can be written as the product 
of two contributions. 



dn{q, gi, g2, q) = dn{q, gi, q) [dn{q, g2, gi) + dn{gi, g2, q) - small correction] (10) 

Emission of two gluons could therefore be described as the emission of one first gluon 
from the original quark anti-quark pair, and the subsequent emission of the second gluon 
of smaller transverse momentum, from one or the other colour dipole involving the doubly 
(colour) charged gluon. In other words, the emission of the first gluon splits the original 
colour dipole into two dipoles which radiate independently. 

The available rapidity range for the emission of the second gluon, which is just the 
sum of the rapidity ranges available in the two dipoles available after the first emission, 
is more than the rapidity range in the original quark anti-quark dipole. 



The quantity Sqg^Sg^q/s is the invariantly defined transverse momentum of the first 
gluon. Eq.(|lT]) therefore implies that the emission of the first gluon increases the available 



phase space for further gluon emission by an amount /n(/c5.i/^i)- 

The above dipole splitting scenario is generalised in the dipole cascade model. Emis- 
sion of the second gluon is thought to give rise to three dipoles which radiate indepen- 
dently, and so on. The available rapidity range for emission of new gluons at a certain 
transverse momentum scale k± is, as a consequence of Eq.(0), a function of all the emis- 
sions above that scale. 



n-l 



^y{gn) = H^) + Y.H^) 



(12) 



The increase of phase space due to gluon emission is conveniently visuahzed using the 
Lund triangular phase space diagrams, such as those shown in Fig.^ The y-axis represents 
the logarithm of the transverse momentum scale k^, denoted by k. Eq.(|^) then means that 
as K is decreased the available phase space increases linearly with —k. This is represented 
by the base line of the triangle which increases as kappa decreases. In the second figure, 
we show one gluon emitted at a certain scale Ki. The phase space for emission at a lower 
scale is the base line of the original triangle plus another term which increases linearly with 
Ki — K according to Eg. (p!TD . This is represented by the additional triangular appendage 
starting at ki. The total phase space is therefore obtained by measuring the length of 
the base line of the structure, going around all triangular appendages. The third figure 
in Fig.|] shows this after several emissions. The forth figure is a view from the "bottom" 
of the phase space triangle. The triangular appendages appear as lines coming out of the 
original base line, increasing its length. 

The transverse momentum scale k is a resolution scale on the partonic state. The 
smaller the scale, the softer the gluons we resolve. Fig.^ illustrates that the available 
phase space scales such that it increases faster than a linear rate when the resolution 
scale is lowered. This is a manifestation of the anomalous dimensions of QCD. 

It should be noted that the phase space in Eq. ( |l2l) is not written in terms of variables 
describing the state from which further emission is being considered. It requires infor- 
mation about the history of the event, i.e. which gluon was emitted first, which came 
next etc., and at what transverse momenta they were emitted. It also ignores the effect 
of recoil due to gluon emissions. For a small invariant mass of the radiating dipole, recoil 
effects could be significant. Emission from any one dipole would give a recoil to the two 
partons in that dipole, and would therefore change the mass of the adjoining dipoles. The 
assumption of independent emission from the dipoles itself, however provides one measure 
defined in terms of the dipoles available for radiation: 

dipoles 

i{kl)= E (13) 

When a string without gluons decays, the rapidity range available to hadrons along 
the mean hyperbolic decay region is given by an expression similar to Eq.(P), with the 
transverse momentum scale kj_, replaced by a suitable scale ml related to hadronisation. 
But generalisation of the phase space available for hadron formation goes along slightly 
different lines. One observes that the hadronisation scale ml is independent of the scale 
at which the partonic state is calculated. Since string fragmentation is infrared stable, 
the phase space for hadrons should be an infrared stable generalisation of the fiat string 
phase space, as softer and softer gluons are added. Addition of gluons softer than the 
scale ml should not affect the phase space significantly. 

Such a generalisation was developed in p, and has been discussed over the years. This 
measure, called the A-measure, is defined in terms of differential equations for a general 
directrix A{^) as follows. We define a four- vector valued functional qr and a scalar valued 
functional T : 

r(0 = i+f;?:^M^.(6-6)T(6) (14) 

Jo mi 
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T(0 Jo 

which satisfy the differential equations: 

dT = 

aqr = dA 5— gr 



qT{0 = T^ l^dA{C)T{a (15) 







dq^ = 2{l -%)qTdA (16) 



ml' 



They can be integrated to: 



qT{i')dA{i') 



no = exp 

qUO = mill -no-'] (17) 



/o mo 



Finally we note that we may define a four-vector valued function = Afj, — q^ so 
that: 

dX = 5 — qr (18) 

Thus the vector g^- is the tangent to the ^Y-curve reaching out to the directrix. From 
Eq.(|l7|) we conclude that the functional T is the exponential of an area (note that the 
vector dA is light-like, and qdA is an area element between the directrix and the A'-curve). 
This area corresponds to the measure A mentioned earlier. We also note that the vector 
qT quickly reaches the length mo. Integrating the differential equations for a directrix 
built up from a set of light-like vectors {kj} we obtain the iterative equations: 

T- 



I3+1 

QTj+i = Ij+iqxj + ^ ' '^^' kj+i (19) 



1 + Ij+i 



The starting values are Tq = 1 and gro = ^o- 

The factorisability of the functional T means that its logarithm, A, can be written as 
a sum containing one term for each kj-. 

n+l n+1 

X = Y.AX, = Y^ In (1 + qT,-ik,/ml) (20) 
i=i i=i 

The AXj defined in this way corresponds to the size of a subarea just as the total A 
corresponds to the area spanned between the X and the directrix curves. In Fig. ^ we 
exhibit such a region or "plaquette" bordered by the "initial" qxj and the "final" qxj+i 
together with a hyperbolic segment from the A'-curve and the gluon energy-momentum 
vector kj from the directrix. 
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Figure 7: A plaquette bordered by qrj, QTj+i , kj and a hyperbolic segment from the 
A'-curve, as described in the text. 



Eq. (|20|) should be compared with Eq. (|T3D . Eq. ( ^) expresses the gluon emission phase 
space as a sum of contributions, one term for each colour dipole. Eq. ( pO[ ) expresses the in- 
frared stable hadronisation phase space as a sum of contributions, one term for each vector 
along the directrix, or one term for each plaquette as mentioned in the previous para- 
graph. With this analogy in mind we called these plaquettes "Generalised Dipoles" (GD) 
in 1^. The hadronisation procedure outlined in Sec. ^ could be thought of as a procedure 
which divides these GDs into smaller plaquettes, one for each hadron. In we discussed 
further the similarities and differences between the A-measure, and the ^-measure. 

In we have also reported on some interesting scaling properties of the contribu- 
tions to the ^-measure from individual dipoles. It was found that for partonic states 
obtained from parton showers in Monte Carlos like Pythia and Ariadne, the dipole 
sizes measured relative to the transverse momentum cut-off scale k± had a distribution 
which was surprisingly insensitive to the scale k± cf. Fig.^. We showed that it is pos- 
sible to understand the distributions in terms of a "gain-loss" analysis: as we go down 
in transverse momentum, the dipole sizes increase, since they are measured relative to 
transverse momentum. But emission of softer gluons would split some dipoles of a given 
size into smaller ones which are lost from a certain bin in dipole size. The same process 
would split some dipoles of larger size, so that they end up in the bin under consideration. 
An integro differential equation taking into account these considerations reproduces the 
observed behaviour in dipole sizes quite well, as illustrated in Fig.||. 

It was also found that the contributions to the A-measure from individual GDs have 
very similar scaling properties, if the scale ml is chosen to be proportional to the ordering 
variable in the parton cascade. Even though there is no direct way to calculate the 
distribution of sizes for GDs, in the light of certain similarities between the A-measure 
and the ^-measure, such behaviour is not completely unexpected. 

That parton showers tend to produce structures of the same size relative to the cut-off 
scale is interesting because the A-measure and the Generalised Dipoles are connections 
between perturbative physics and hadronisation. On the one hand they are related, and 
share many properties with phase space for gluon emission. On the other hand the 
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Figure 8: Dipole sizes obtained from Pythia and Ariadne , measured relative to the 
cut-off scales, are relatively insensitive to the cut-off scales. The behaviour is quite well 
described by an MLLA calculation. 



fragmentation process reviewed in Sec.^ is a connection between the plaquettes we called 
GDs and the plaquettes we defined for each hadron in Sec.|[ Properties of the colour 
dipoles which are preserved in the Generalised Dipoles will be carried over to the final 
state hadrons. 

4 Conclusions 

String fragmentation can be formulated as the production of a set of "plaquettes" between 
the partonic curve or directrix, and the hadronic curve. The sum of the areas of the 
plaquettes is taken as the area in the Lund Model area law. The mean size of this total 
area can be calculated from the specified perturbative state, and is proportional to the 
length of a curve called the A'-curve. The length of the A'-curve is also called the A- 
measure. The A-curve looks like a set of connected hyperbolae stretching across planar 
regions around the directrix which we called "Generalised Dipoles" . The size of the GDs 
measured relative to the parton shower cut-off scale seems to have a distribution which 
is insensitive to the cut-off scale itself. This can be understood as a balance between 
emission of more gluons at a softer scale versus the scaling of the phase space variable 
itself. 
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